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^SJ ' Two examples of recent progress in applications of the Dyson-Sch winger equation 

^~^ . (DSE) formalism are presented: 

—1. I (1) Strong coupling quantum electrodynamics in 4 dimensions (QED4) is an often studied 

OS ' model, which is of interest both in its own right and as an abelian model of quantum chro- 

^ . modynamics (QCD). We present results from a study of subtractive renormalization of the 

Qh! fermion propagator Dyson-Schwinger equation (DSE) in massive strong-coupling quench- 

Q_i' ed QED4. Results are compared for three different fermion-photon proper vertex Ansdtze: 

^ • bare 7^^, minimal Ball-Chiu, and Curtis- Pennington. The procedure is straightforward to 

implement and numerically stable. This is the first study in which this technique is used 

.5^ I and it should prove useful in future DSE studies, whenever renormalization is required in 

^ ■ numerical work. 

(2) The Bethe-Salpeter equation (BSE) with a class of non-ladder scattering kernels is 
solved in Minkowski space in terms of the perturbation theory integral representation 
(PTIR). We consider a bound state of two spinless particles with the formal expression of 
the full scattering kernel in a cp'^a scalar model. Making use of the PTIR we isolate the 
possible kinematical singularities of the bound state amplitude in Minkowski space. The 
resulting BS amplitude is written as a parametric integral of its weight function in this 
approach. We derive an integral equation for the weight function with a real kernel. We 
compare numerical solutions of the non-ladder scattering kernel with that of the massive 
scalar exchange kernel in the ladder approximation. 

1 Quenched Massive QED4 

Strong coupling QED in three space and one time dimension has been studied within the 
Dyson-Schwinger Equation (DSE) formalism for some time |jl|, ^, Q, in order to see whether 
there may be a phase transition to a nontrivial "local" theory at high momenta 0, 0, as a 
model for dynamical chiral symmetry breaking (DCSB) in walking technicolor theories |^ |^, 
and also as an abelianized model for nonperturbative phenomena in QCD 0, |^. For a recent 
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review of Dyson- Schwinger equations and their application see Ref. 0. The usual approach is 
to write the DSE for the fermion propagator or self-energy, possibly including equations for the 
photon vacuum polarization [0, |11| or the fermion-photon proper vertex |^. An appropriate 



Ansatz is made for the undefined Green's functions that contain the infinite continuation of the 
tower of DSE's. The resulting nonlinear integral equations are converted to Euclidean metric 



in the usual way [0 and solved numerically by iteration from an initial guess. DCSB occurs 
when the fermion propagator develops a nonzero scalar self-energy in the absence of an explicit 
chiral symmetry breaking (ECSB) fermion mass. We refer to coupling constants strong enough 
to induce DCSB as supercritical and those weaker are called subcritical. We write the fermion 
propagator as 

Z{p') I 

with Z{p^) the finite momentum- dependent fermion renormalization, and Bijp') the scalar self- 
energy. In the massless theory (i.e., in the absence of ECSB) by definition DCSB occurs when 
S(/) ^ 0. Note that A{p'^) = l/Z{p'^) and M(p2) = B (p^) / A{p'^) . 

Many studies, even until quite recently, have used the bare vertex as an Ansatz for the 



S(P) = -u ^JTZ^ = ^uZ^IVZl Y^TZ^ (1) 



one-particle irreducible (1-PI) vertex ^{k^p) [0, ||, ||, |^, |1^, ^ |T3|, despite the fact that this 
violates the Ward-Takahashi Identity (WTI) [l^. It is also common, especially in studies 
motivated by walking technicolor theories [^], to find vertex Ansdtze which claim to solve the 
WTI, but which still possess kinematic singularities in the limit of zero photon momentum 
g^ = {k ~ p)'^ ^ ||15|, |l6l. With any of these Ansdtze the resulting fermion propagator is not 
gauge-covariant, i.e., physical quantities such as the critical coupling for dynamical symmetry 
breaking, or the mass itself, are gauge-dependent |jl5|. A general form for r''(A;,p) which does 



satisfy the Ward Identity was given by Ball and Chiu in 1980 |T^; it consists of a minimal 
longitudinally constrained term which satisfies the WTI, and a set of tensors spanning the 
subspace transverse to the photon momentum q. 

Although the WTI is necessary for gauge-invariance, it is not a sufficient condition; further, 
with many of these vertex Ansdtze the fermion propagator DSE develops overlapping logarithms 
so that it is not multiplicatively renormalizable. There has been much recent research on the 
use of the transverse parts of the vertex to ensure both gauge-covariant and multiplicatively 



renormalizable solutions [1S]^|23|, some of which will be discussed below. 

What is common to all of the studies so far is that the fermion propagator is not really 
subtractively renormalized. Most of these studies have assumed an initially massless theory 
and have renormalized at the ultraviolet cutoff of the integrations, taking Zi = Z2 = 1. Where 
a nonzero bare mass has been used 0, |T^ , it has simply been added to the scalar term in the 



propagator. Although there have been formal discussions of the renormalization 0, |2T|, the 
important step of subtractive renormalization has not been performed. 

1.1 DSE, and Vertex Ansatze 

The DSE for the renormalized fermion propagator, in a general covariant gauge, is 

S-\p')=Z,{^,,A)[^~m'{A)]-tZ,{ii,A)e'f-0^^rS{k)V^{k,p)D,M; (2) 

here q = k — p is the photon momentum, /i is the renormalization point, and A is a regular- 
izing parameter (taken here to be an ultraviolet momentum cutoff). We write mo(^) for the 



regularization-parameter dependent bare mass. In the massless theory (i.e., in the absence of 
an ECSB) the bare mass is zero, mo (A) = 0. The physical charge is e (as opposed to the bare 
charge eo), and the general form for the photon propagator is 



with ^ the covariant gauge parameter. Since we will work in the quenched approximation and 

--0 



the Landau gauge we have e^ = Cq = 47rao and 



D'^^(g)-D(;^(g)=(-^'^^ + ^U, (4) 



for the photon propagator. 

1.1.1 Vertex Ansatz 

The requirement of gauge invariance in QED leads to the Ward-Takahashi Identities (WTI); 
the WTI for the fermion-photon vertex is 

q,V^{k,p) = S-\k)-S-\p), (5) 

where q = k — p. This is a generalization of the original differential Ward identity, which 
expresses the effect of inserting a zero-momentum photon vertex into the fermion propagator, 

—- — = r {p,p) . 6 

OPu 

In particular, it guarantees the equality of the propagator and vertex renormalization con- 
stants, Z2 = Zi. The Ward-Takahashi Identity is easily shown to be satisfied order- by-order in 
perturbation theory and can also be derived nonperturbatively. 

As discussed in 0, this can be thought of as just one of a set of six general requirements 
on the vertex: the vertex must satisfy the WTI; it should contain no kinematic singularities; it 
should transform under charge conjugation (C), parity inversion (P), and time reversal (T) in 
the same way as the bare vertex, e.g., 

C^T^{k,p)C = -Tl (7) 

(where the superscript T indicates the transpose); it should reduce to the bare vertex in the 
weak-coupling limit; it should ensure multiplicative renormalizability of the DSE in Eq. (^; 
the transverse part of the vertex should be specified to ensure gauge-covariance of the DSE. 
Ball and Chiu |]T^ have given a description of the most general fermion-photon vertex that 



satisfies the WTI; it consists of a longitudinally-constrained (i.e., "Ball-Chiu") part EgQ, which 
is a minimal solution of the WTI, and a basis set of eight transverse vectors T^{k^ p) , which span 
the hyperplane specified by q^Tj^{k,p) = 0, q = k — p. The minimal longitudinally constrained 
part of the vertex is given by 

rUk,p) = llAik^) + A{p^)]Y + ^^ l[A{k^) - A{p^)]^ - [B{k^) - B{p^)]\ . (8) 



A general vertex is then written as 

T'^ik,p) = V^scik^p) + J:n{k\p\q')Tt{k,p) , (9) 

where the Tj are functions which must be chosen to give the correct C, P, and T invariance 
properties. 



The work of Curtis and Pennington |^ was mentioned above in connection with the speci- 
fication of the transverse vertex terms in order to produce gauge-invariant and multiphcatively 
renormahzable solutions to the DSE. In the framework of massless QED4, they eliminate four of 
the transverse vectors since they are Dirac-even and must generate a scalar term. By requiring 
that the vertex r^(/c,p) reduce to the leading log result for k ^ p they are led to eliminate 
all the transverse basis vectors except Tq, with a dynamic coefficient chosen to make the DSE 
multiphcatively renormahzable. This coefficient has the form 

T,ie,p\q') = ^[Aik')-Aip')]/dik,p) , (10) 

where d{k,p) is a symmetric, singularity free function of k and p, with the limiting behavior 
Iimfc2^p2 (i(fc,p) = k"^. [Here, v4(p^) = 1/Z(p^) is their l/jF(p^).] For purely massless QED, 
they find a suitable form, d{k,p) = [k"^ — p^)^/(fc^ +P^)- This is generalized to the case with a 
dynamical mass M{p'^), to give 

They establish that multiplicative renormalizability is retained up to next-to-leading-log order 
in the DCSB case. Subsequent papers establish the form of the solutions for the renormalization 
and the mass |]2T[, and demonstrate the gauge- covariance of the solutions |^. (However, Dong, 



Munczek and Roberts [22] have recently contested this). Bashir and Pennington |2^ have 



recently described a vertex Ansatz which makes the fermion self-energy exactly gauge-covariant, 
in the sense that the critical point for the chiral phase transition is independent of gauge. 

In this work we will primarily compare the Curtis- Pennington Ansatz with results using the 
bare vertex. Some solutions are also obtained with the minimal Ball-Chiu vertex which, like 
the Curtis-Pennington vertex, satisfies the WTI, but which does not lead to gauge-invariant 
solutions. 

For each vertex Ansatz, the equations are separated into a Dirac-odd part describing the 
finite propagator renormalization A{p'^), and a Dirac-even part for the scalar self-energy, by 
taking ^Tr of the DSE multiplied by ^/p^ and 1, respectively. These equations are rotated to 
Euclidean metric, giving equations for the spacelike momenta only. The volume integrals / d'^k 
are separated into angle integrals and an integral / dk"^; the angle integrals are easy to perform 
analytically, yielding the two equations which which will be solved numerically. 

1.2 The Subtractive Renormalization 

The subtractive renormalization of the fermion propagator DSE proceeds similarly to the one- 
loop renormalization of the propagator in QED. (This is discussed in |p and standard texts.) 
One first determines a finite, regularized self-energy, which depends on both a regularization 



parameter and the renormalization point; then one performs a subtraction at the renormal- 
ization point, in order to define the renormahzation parameters Zi, Z2, ^3 which give the full 
(renormalized) theory in terms of the regularized calculation. 

A review of the literature of DSEs in QED shows, however, that this step is never actually 
performed. Curtis and Pennington |2I[] for example, define their renormalization point at the 



UV cutoff. Miransky [Q gives a formal discussion of the variation of the mass renormalization 
Zm(/i, A), but does not implement it numerically. 

Many studies take Zi = Z2 = 1 [0, |20|, ^; this is a reasonable approximation in cases 
where the coupling oq is small (i.e., Oq ~ 1-15), but if Oq is chosen large enough, the value of 
the dynamical mass at the renormalization point may be significantly large compared with its 



maximum in the infrared. For instance, in [^, figures for the fermion mass are given with Oq 



= 0.97, 1.00, 1.15 and 2.00 in various gauges. For ao = 2.00, the fermion mass at the cutoff is 
down by only an order of magnitude from its limiting value in the infrared. 

Repeating the arguments from [0, one defines a regularized self-energy S'(/i, A;p), leading 
to the DSE for the renormalized fermion propagator, 

S-\p) = Z2(^,A)[^-mo(A)]-S'(/x,A;p) 

= ^-m(^) -S(^;p) , (12) 

where the (regularized) self-energy is 

S'(/i, A;p) = zZi(/x, A)e' |^ -^T^Sifi; k)f^ifi; k,p)D,,ifi; (p - k)) . (13) 

[To avoid confusion we will follow Ref. [^] and in this section only we will denote regularized 
quantities with a prime and renormalized ones with a tilde, e.g. S'(yU, A;p) is the regularized 
self-energy depending on both the renormalization point /i and regularization parameter A 
and S(/i;p) is the renormalized self-energy.] As suggested by the notation (i.e., the omission 
of the A-dependence) renormalized quantities must become independent of the regularization- 
parameter as the regularization is removed (i.e., as A ^ 00). The self-energies are decomposed 
into Dirac and scalar parts, 

S'(/i, A; p) = E:,(/i, A; p2) ^ + E',(/i, A; p^) 

(and similarly for the renormalized quantity, S(/x,p)). By imposing the renormalization bound- 
ary condition, 

S-\p) , =^-m(/i), (14) 



one gets the relations 

for the self-energy, 

for the renormalization, and 



S,,,(/^; p') = S',_,(/i, A; p') - S:,^,(/i, A; /x^) (15) 

Z2(/i,A) = l + S:^(/i,A;/i2) (16) 



mo(A) = [m{i2) - S',(^, A; ^^)J /Z2(/x, A) (17) 

for the bare mass. In order to reproduce the case with no ECSB mass, for a given cutoff A, one 
chooses m{fi) = S'j,(/i, A;/x^) so that the bare mass rTio(A) is zero. The mass renormalization 
constant is given by 

Z^(/i,A)=mo(A)/m(/i), (18) 



i.e., as the ratio of the bare to renormahzed mass. 

The vertex renormahzation, Zi(/x, A) is identical to Z2(/x, A) as long as the vertex Ansatz 
satisfies the Ward Identity; this is how it is recovered for multiplication into S'(yU, A;p) in 



Eq. (13). It will be noticed that this is inappropriate for the bare-vertex Ansatz since it fails 
to satisfy the WTI; nonetheless, since for the bare vertex case there is no way to determine 
Zi(yU, A) independently we will use Zi = Z^ for the sake of comparison. In the Landau gauge 
for the bare vertex these will then both be 1, since in this case S^(/x, A;p^) = for all p^ as is 
well known 0. 

1.3 Results 

Solutions were obtained for the DSE with the Curtis-Pennington and bare vertices, for cou- 
plings ao from 0.1 to 1.75; solutions were also obtained for the minimal Ball-Chiu vertex, with 
couplings ao from 0.1 to 0.6 (for larger couplings the DSE with this vertex was susceptible 
to numerical noise). In Landau gauge, the critical coupling for the DSE with bare vertex is 
^bare _ ^^g. i^^ critical coupliug for the Curtis-Pennington vertex is a^ = 0.933667 |^, and 



that for the Ball-Chiu vertex is expected to be close to these two values. A full discussion of 
the numerical results with detailed figures can be found elsewhere [^. A brief summary of 
these is given here in the conclusions. 

2 Minkowski- space Bethe-Salpeter Equation 

Considerable interest has been recently attached to the covariant description of bound states in 
conjunction with model calculations of high-energy processes, such as deep inelastic scattering. 
A fully covariant description of composite bound states is essential. The bound state nature 
of hadrons is described by the appropriate vertex function, or equivalently the Bethe-Salpeter 
(BS) amplitude. In a relativistic field theory two-body bound states are described by the 



Bethe-Salpeter (BS) amplitude |2§]. The BS amplitude obeys an integral equation whose kernel 
has singularities due to the Minkowski metric. The resultant solutions are mathematically no 
longer functions but "distributions". This singular structure makes it difficult to handle the 
BS equation numerically. In order to handle such a singular integral equation the analytic 
continuation of the relative-energy variable, which is called "Wick rotation", is widely used 
T7| . The ladder BS amplitude is solved as a function of Euclidean relative momentum in 
the standard approach. If one uses a "dressed" propagator for constituent particles or more 
complicated kernels in the BS equation, the validity of the Wick rotation becomes highly 
nontrivial, e.g., almost all the dressed propagator studied previously in the Dyson-Schwinger 
equation approach contains pathological complex "ghost" poles[^. It is therefore preferable to 
formulate and solve the BS equation directly in Minkowski space. 

We present a method to solve the BS equation without Wick rotation by making use of 
the perturbation theory integral representation (PTIR) for the BS amplitude ||28||. This integral 



representation has been studied for a scalar-scalar bound state in the ladder approximation p9[| . 
We extend this method to a wide class of non-ladder kernels. We rewrite the BS equation as 
the integral equation of the weight function for the BS amplitude and discuss the singularity 
structure of the kernel function for the weight function. 



2.1 Scalar- Scalar BS Equation 

Let us consider a bound state of two spinless particles (j) having a mass m. They interact each 
other through the exchange of another spinless particle a with a mass fi. Let the interaction 
between (p and a be the Yukawa couphng: £j„t = —gcpa. The Bethe-Salpeter amplitude 
$(p, P) for the bound state having the total momentum P and the relative one 2p obeys the 
following equation: 

[D{p/2 + p)D{p/2 - p)]-'^p, P) = j j^Ap, q; PMq, p) (19) 

where D{q) is the propagator of ^-particle and we approximate it with the tree one: D^{q) = 
l/(m^ — q^ — ie). The scattering kernel /(p, g; P) describes the process: 0i02 — >■ 0304 and the 
momentum 2p and 2g are the relative momentum of initial and final states. We consider the 
following formal expression for the full scattering kernel: 



I{p,q;P) = jd^jdi ^ 



le 



A ^ 

+ (cychc permutation of s, t, u) (20) 

where qi is the momentum carried by (^i and s,t and u are Mandelstam variables. The 
symbol A denotes the integral region of 6 dimensionless Feynman parameters ^j such that 
A = \Xi\ii > 0, Yl,ii = 1 (^ = 1) ■ ■ • ) 6)}. The "mass" parameter 7 represents a spectrum of 
the scattering kernel. The function Pch{l,^',g) gives the weight of the spectrum in a different 
channel; ch = {st}, {tu}, {us}. This expression has been derived by Nakanishi and is called 
the perturbation theory integral representation (PTIR)|^^. It should be mentioned that any 
perturbative Feynman diagram for the scattering kernel can be written in this form, so that 
the weight function pi is, in principle, caluculable as a power series of coupling constant g. 

2.2 PTIR for BS amplitude 

Let us consider the s-wave bound state for simplicity^ We assume that the BS amplitude 
<l>(p, P) has an integral representation of the form: 

00 1 



$(p,P) = -i da dz ^"^ ' ' —2 (21) 



-co 



where the non-negative integer parameter ra is a dummy parameter, since a partial integration 
with respect to a changes the power. We can utilize this artificial degree of freedom for a 
numerical study. Substituting the above expression into Eq.(|T^ together with Eq. (ppj) , we 
obtain the following integral equation for </?„(«, z) as follows: 

00 1 00 ( ^. \ 

(pn{a, z) = / da dzJ2 d-f d^ "^(4' )2 — ^ri{a, z; a, z) ipn{oi, z), (22) 



ch 



^ Extension to higher partial wave solutions is straightforward. 



where we have suppressed the dependence on the kernel parameters. Note that the Eq.p^) is 
frame-independent. The real kernel function Kn{(y,z;a,z) with a fixed kernel parameter set 
(7, ^) has the following structure: 

d _ _ 

Kn{a,z;a,z) = — {a^9{a)) hn{a, z) — 9 {{a — uji{a,z,z)){a — uj2{a,z, z))) 

^ J gn{a,z;a,z) , 7 /- - nI /oqn 

x< _ _ _ ^ + kn{a,z]a,z)> (23) 

[yj{a-ui{a,z,z)){a-uj2{a,z,z)) J 

where uji{a, z, z), gn{o:, z; a, z) and /in(«, z) are regular functions. The function kn{a, z\ a, z) is 
also regular in the simple one-a-exchange ladder kernel, but in general it contains a singularity 
such as Pf ■ l/(a — r(a,^,2;))" where the symbol Pf- stands for Hadamard's finite part and 
ripL^ z, 2;) is a regular function. Thus the first term of Eg . (^3]) contains a (5-function singularity 
at a = 0, only if n = 0. This singularity, independent of z^a and z, corresponds to the pole 
singularity in $(p, P) which comes from the free propagation of two (/)'s. On the other hand, 
the second term contains a discontinuity due to the step function depending on a, z, a and 
z. In addition to the singularity in kn{pc^ z\ a, z) as mentioned above, the term has a square 
root line singularity at the boundary of its support. Since the Hadamard's finite part Pf ■ l/x 
coincides Cauchy's principal value, the kernel function i^„=i(Q;, z; a, z) is integrable, so that the 
integral equation (^) with a constant kernel parameter set is numerically tractable by setting 
the dummy parameter n = 1, provided that the weight function Lpi{a,z) is differentiable at 
the singular points of Ki{a^ z\ a, z). As a check on our method we have reproduced the known 
results for the one-a-exchange ladder kernel by solving Eq.(p2D numerically. The numerical 
solution has been extended also to a class of non-ladder ("generalized ladder") kernels. The 
case for the most general scattering kernel, whose weight function can contain a derivative of 
5-function, is currently under investigation. 

3 Conclusions 

We have described preliminary results in a study of four-dimensional quenched QED, with sub- 
tractive renormalization performed numerically, on-the-fly during the calculation. We believe 
that this is the first calculation of its kind, and the technique described here will be applicable 
elsewhere (e.g., in both QED and QCD), whenever numerical renormalization is needed. 

The Curtis-Pennington vertex has been the primary focus of this study, since it has the 
desirable properties of making the solutions approximately gauge-invariant and also multi- 
plicatively renormalizable up to next-to-leading log order. Solutions have been obtained for 
comparison purposes, using the minimal Ball-Chiu vertex and using the bare vertex Ansdtze 
(with Zi = Z2). For the Ball-Chiu vertex, couplings in the range from a^ = 0.1 to 0.6 were 
used, while couplings up to ao = 1-75 were used with both the bare and the Curtis-Pennington 
vertices. Various renormalization points and renormalized masses were studied. 

The subtractive renormalization procedure is straightforward to implement. The solutions 
are stable and the renormalized quantities become independent of regularization as the reg- 
ularization is removed, which is as expected. For example, the mass function M(p^) and 
momentum-dependent renormalization ZijP') are unchanged to within the numerical accu- 
racy of the computation as the integration cutoff is increased by many orders of magni- 
tude. For the range of couplings ao considered, the values of the renormalization constant 
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Zi{fi,A) = Z2{fi,A) are never very far from 1 and vary relatively weakly with the choice of 
renormahzed mass and cutoff as expected in Landau gauge. For subcritical couplings, us- 
ing the Curtis-Pennington vertex, we also find that the mass renormalization Zm(/^, A) scales 
approximately as Z^(yU, A) oc (/iVA^)^5-7cp(ao))^ where e.g. 7cp(0.5) = 0.358. 

An unexpected feature of the equations is that for any set {ao, /i^, 'm(/i)}, as the cutoff is 
increased there is a region where the dynamical mass is negative. The significance of this for 
QED is not completely understood. One possibility is that it may signal the the failure of 
multiplicative renormalizability for the model DSE, in which case further refinements of the 
vertex Ansatz may be called for. It should be emphasized that this negative dip in the scalar 
self-energy is very small unless the coupling oq approaches 2. For instance, for ao = 1.00, the 
negative peak was ~ 2.4 x 10~^ the size of the renormahzed mass. 

Extensions of this work to include other gauges and vertices of the Bashir-Pennington type 
are underway [|^ . 

For the Bethe-Salpeter equation we have derived a real integral equation of the weight 
function for the scalar-scalar BS amplitude with a formal expression of the full scattering 
kernel. We found that our integral equation is numerically tractable for a class of non-ladder 
scattering kernels. We have verified that our numerical solutions agree with those previously 
obtained from a Euclidean treatment of the pure ladder limit. This represents a powerful new 
approach to obtaining solutions of the BSE and a more detailed discussion will appear soon[^ . 
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